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We discuss transport and localization properties on the insulating side of the disorder dominated superconductor- 
insulator transition, described in terms of the dirty boson model. Analyzing the spectral properties of the 
interacting bosons in the absence of phonons, we argue that the Bose glass phase admits three distinct 
regimes. For strongest disorder the boson system is a fully localized, perfect insulator at any temperature. 
At smaller disorder, only the low temperature phase exhibits perfect insulation while delocalization takes 
place above a finite temperature. We argue that a third phase must intervene between these perfect insula- 
tors and the superconductor. This conducting Bose glass phase is characterized by a mobility edge in the 
many body spectrum, located at finite energy above the ground state. In this insulating regime purely elec- 
tronically activated transport occurs, with a conductivity following an Arrhenius law at asymptotically low 
temperatures, while a tendency to superactivation is predicted at higher T. These predictions are in good 
agreement with recent transport experiments in highly disordered films of superconducting materials. 
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1 Introduction 

The disordered superconductor-to-insulator transition (SIT) HI |2l has recently attracted a lot of attention 
due to a series of surprising and puzzling transport experiments, especially on the insulating side of the 
quantum phase transition. In strongly disordered samples with nearly critical disorder, a large magne- 
toresistance peak is observed as the sample is tuned further away from the superconducting phase by a 
magnetic field IS). The latter HUH and several further experimental indications in transport |l5l|6l and 
tunneling 111 are commonly interpreted as an indication for the survival of an electronic pairing gap close 
enough to the transition, despite the loss of global phase coherence in this kind of samples. This scenario is 
supported by numerical studies which suggest that the pairing gap and the existence of bosons (pairs) can 
survive both the disorder driven lH and the field-driven SIT ||9l (in the absence of Coulomb interactions). 
However, apart from the above described type of SI transitions, there is another class of systems with 
weaker disorder ||2l I?) in which superconductivity is predominantly suppressed by Coulomb interactions 
and the concomitant destruction of electron pairing. 

The experimental observations in strongly disordered samples motivate the assumption that the low 
energy physics in the vicinity of the SIT can be described by interacting bosons subject to a disorder 
potential. This forms a cornerstone of the dirty boson model for the SIT 1 10||, which predicts that disorder 
(or a magnetic field) drive the superconducting state into an insulating, but compressible Bose glass phase. 
In sufficiently pure samples with a regular granular structure and commensurate filling a transition to a 
Mott insulating phase with a hard gap for charged excitations is predicted to exist at smaller hopping 
strength. However, such an incompressible phase has little reason to occur in strongly disordered thin 
films, independent of the presence or absence of self-generated granularity. 

In this paper we adopt a purely bosonic description of the phase transition for dimensions d > 2. We 
focus on the properties of intrinsically electronic transport in the insulating Bose glass phase, where we 
strictly neglect electron-phonon coupling. The latter is empirically known to be rather weak in the relevant 
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experimental systems iTTl. The behavior of the intrinsic conductivity turns out to be surprisingly rich. Our 
analysis leads to a generic picture of the physics close to the SIT which we believe to be key to understand- 
ing several unexplained features of transport in Bose glasses. On one hand, it has been a long standing 
puzzle that strongly disordered superconducting films exhibit a simply activated (Arrhenius) resistance at 
low temperatures lfT2l instead of variable range hopping, which is normally expected in disordered insu- 
lators. While the absence of variable range hopping may be blamed on the weak electron-phon coupling, 
most simple explanations for an Arrhenius resistance must be ruled out: The presence of a hard gap for 
bosons ifTSl is excluded in the strongly disordered, compressible Bose glass. Further, an explanation based 
on nearest neighbor hopping is inconsistent with the observed magnitude of the pre-exponential factor of 
the conductivity. Finally, interpreting the activation energy as the energy to form unpaired electrons lfT4l is 
ruled out by the observed positive magnetoresistance. 

Another interesting transport feature was discovered very recently. At ultralow temperatures, Bose 
glasses such as InO^ [[3] or TiN fT5\ exhibit huge resistance jumps when driven beyond a certain threshold 
in voltage [TT|. This phenomenon was explained as a thermal run-away into an overheated state of the 
electronic subsystem which is nearly decoupled from the phonons lfT6l . This interpretation relies crucially 
on the assumption that the resistance is sensitive to the electron temperature rather than the lattice temper- 
ature. As stressed in Ref. I1I6I . this is at odds with the standard mechanism of phonon assisted transport 
in insulators and implies that the transport in the Bose glass phase must essentially be activated by the 
electrons themselves. However, no microscopic mechanism was proposed. 

In this paper, we argue that purely electronic transport and its simply activated form are both natural to 
expect close enough to the SIT for d > 2, provided that electron-phonon coupling plays a negligible role. 

2 The Bose glass 

For the general analysis of the phase diagram below, a specific microscopic picture of disordered bosons 
is not really essential. Nevertheless, in order to fix ideas it is helpful to argue in terms of the model 
considered by Ma and Lee ifTTl . Let ipi be the Anderson localized states of non-interacting electrons, each 
with a localization length of typical size ^. We then restrict the Hilbert space of the full Hamiltonian (with 
an attractive BCS interaction) to configurations with all states i either empty or doubly occupied with a 
spin singlet. This reduces the problem to a system of disordered hard core or pseudo-spin- 1/2 bosons ifTSll : 

H = - Y^ie. - ^i)sl + ^(tys+^j + l^-^-)- (1) 

Each state i has a random energy e;, the average density of states being 1^. The level spacing in the 
localization volume of bosons is thus 6^ = 1/ Sf^v, and ^ denotes the chemical potential. In the Hamiltonian 
([T]) we have omitted a density-density interaction j -^ij^i^j' would arise, e.g., from Coulomb 
interactions. However, it would not affect our general considerations as long as Jij are short ranged and 
weak enough so as not to introduce frustration and glassy effects. 

Note that the attraction-induced hopping of bosons, ty , is appreciable only between spatially overlap- 
ping states, tij cxp[— jr; — |/^], with the center of mass coordinates rij of the states. If the states i, j 
are within the same localization volume, ty will typically decrease as a power law with distance in energy 
ei — ej. The superconductivity in such a model, built from nearly critical wavefunctions tpi, has recently 
been studied in great detail fT9l. 

Let t be the typical magnitude of tij corresponding to nearest neighbors in space and energy. The 
dimensionless parameter g = characterizes the strength of disorder as compared to the kinetic energy. 
If g ^ 1, the disorder plays very little role and the physics will be equivalent to that of an infinite-?/ 
Hubbard model at incommensurate filling: The bosons delocalize and condense into a superconducting 
state with infinite conductivity below a critical temperature Tc{g) which decreases with growing g, see 
Fig.[T] However, at some critical gc = 0{1) the disorder becomes dominant and localizes the bosons, 
transforming the superconductor into an insulating Bose glass. The associated quantum phase transition is 
expected to be of second order ifTOll . 
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Fig. 1 Phase diagram of the SIT and localization in the Bose glass. At strong enough disorder, g > g^, 
superconductivity is suppressed and gives way to the insulating Bose glass. The T = Q mobility edge, 
Ec, rises monotonously from zero at criticality and leads to An'henius behavior of the conductivity at very 
low temperature. This part of the phase diagram should be compared to the data of 1121 . At stronger 
disorder, g > g,, non-extensive excitations localize at T = 0, the collectively localized phase with a = 
being stable up to a finite Tioc where a phase transition to a conducting phase occurs. At strongest disorder 
g > goo the localized phase persists up to infinite temperature. 

2. 1 Spectrum of the Bose glass 



The transport properties of an interacting insulator can best be understood by analyzing its spectrum 1201 
|2Tl|22l. If the excitation spectrum is discrete (point spectrum) at low energies, the system is localized 
and diffusion is absent. Indeed, the discreteness implies that a local excitation from an initial stationary 
state only couples to a finite number of eigenstates of the full system. The subsequent dynamics is thus 
quasiperiodic and non-diffusive, the excitation remaining local. In such a state thermalization cannot take 
place after the excitation, the system being in a genuine quantum glassy state ||2T]| . In contrast, a continuous 
spectrum reflects the fact that local excitations couple to infinitely many eigenstates of the system, so that 
diffusion, conduction and thermalization are possible. 

It is clear that on the superconducting side of the transition the spectrum is a continuum. Indeed, the 
spontaneous breaking of the continuous U{1) symmetry and the associated long range order ensure the 
presence of delocalized Goldstone modes down to the lowest energies in dimensions d > 2. Their delocal- 
ized nature implies a continuum in the spectrum of excitations at all temperatures. On the Bose glass side 
of the transition, the continuous spectrum cannot persist down to w = 0, at least at T = 0. Indeed, if there 
were delocalized modes at arbitrarily low energies, one would expect Bose condensation into these modes, 
in contradiction to the absence of superfluidity. If condensation were avoided despite delocalization, one 
would have a Bose metal, a phase commonly believed not to exist in short range interacting bosons. 1231 
Hence, if extended low energy states can be excluded, the phase beyond gc must be an insulating Bose 
glass with a discrete spectrum and localized excitations at energies below some threshold Ec{g). On the 
other hand, assuming the SIT to be continuous, the mobility edge Ec separating the discrete spectrum from 
the continuum must evolve continuously with g. In particular, it must vanish at criticality, Ec {g i .gc) — > 0. 

It is important to establish whether or not Ec is an extensive energy l2TI . We argue that close enough 
to the SIT the mobility edge Ec is finite, while it diverges and becomes extensive at some larger disorder 

> 5c- Consider the addition of an extra boson to the ground state of the Bose glass, whereby the boson 
is injected at an energy uj above the chemical potential ^. At small enough lo, the extra boson cannot 
delocalize because there are no extended states at the chemical potential in the Bose glass. However, one 
expects that if w ;2> t, the extra boson creates an excitation which does delocalize, provided g — gc 'is small 
enough. Indeed, at 5 = gc we know that bosons at the chemical potential delocalize and establish long 
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range phase coherence. This happens despite their mutual repulsion and the reduced phase space available 
to each boson. On the other hand, the phase space to hop and delocalize is significantly larger for a boson 
at higher energies uj ^ t. Collisions with other bosons are indeed rare beyond the strip of energies within 
/i ± i where quantum fluctuations in the ground state are significant. Therefore, those high energy states 
will remain extended beyond the disorder strength g — gc, where only the more restricted states at w = 
start localizing. Clearly, the excitation energy of these particular delocalized states are an upper bound for 
the mobility edge Ec, implying that Ec remains non-extensive close enough to the SIT. 

The above reasoning applies strictly only for d > 2 where single bosons always delocalize in weak 
enough disorder, or to d = 2 with spin-orbit scattering potentials for spinful bosons. One may worry that 
the situation is different in = 2 without spin-orbit scattering, since single particles always localize (see 
Ref. 1281 ). However, localization in weak disorder in 2d is due to very subtle interference effects, which 
are likely to be destroyed by the remaining interaction effects with other bosons, even if the injection of an 
extra boson a high energy is considered. We expect this to preempt the scenario of Ec jumping immediately 
to extensive values at gc in d — 2, which would imply the simultaneous localization of all extended states 
at intensive energy. 

If the disorder strength is increased beyond bosons injected at arbitrarily high energy will become 
localized, as well as any local excitation of non-extensive energy. This is indeed expected using arguments 
similar to Refs. 12411211 . Once the level spacing for bosons becomes significantly larger than their inelastic 
scattering rate from other bosons, a localization transition in Fock space is predicted. A lower bound 
for can be obtained from the criterion that all single boson injections be localized at arbitrarily high 
energies cj ^ t above /i. Since interaction effects are weak at these high energies, g* can be estimated 
by the location of complete Anderson localization of single bosons which float energetically far above the 
chemical potential of the ground state. Such an estimate provides a (presumably rather tight) bound for g^. 

2.2 Finite temperature 

The above arguments only allowed us to conclude that the conductivity at T = vanishes in linear re- 
sponse, as it happens by definition in any insulator. The difference between the regimes of disorder dis- 
cussed above is much more dramatic at finite temperature. If the lower edge Ec of the continuum is finite, 
there is a finite density ~ cxp(— E'c/T) of excitations above Ec at asymptotically small temperatures. 
Low energy excitations which were localized at T = inelastically collide with those excitations in the 
continuum and thereby acquire an inverse life time (level width). The latter is proportional to the collision 
probability and thus essentially to cxp{—Ec/T). The whole spectrum thus becomes continuous immedi- 
ately, see Fig. |2l and the conductivity will be finite for all T > 0. Under these conditions, the simplest 
conduction channel consists in the diffusion of thermally activated excitations above the mobility edge, 
provided the latter carry charge. Alternatively, it is also possible to have hopping transport of low energy 
bosons, whereby energy conservation is assured by the emission or absorption of modes above the mobil- 
ity edge. At very low temperature both channels have a bottle neck arising from the exponentially small 
probability of finding excitations above Ec- Accordingly, the conductivity is expected to be asymptotically 
of AiThenius form. More precisely, if we write 

a(r, g) = ao exp[-^^(.g, T)/T] (2) 

one expects the effective activation energy to tend to ^^^(g < gc,T ^ 0) ^ Ec{g)- Note that the above ar- 
guments assume that the mobility of excitations does not strongly decrease as E I Ec, and that the spectral 
weight in the continuum is not too inhomogeneously distributed in energy and space. If neither assumption 
holds, there can be significant, but subleading, temperature corrections to the Arrhenius behavior. 

Note that this conduction mechanism is entirely electronic in nature and does not depend on phonons. 
It thus provides a possible rationale for the occurrence of thermal run-away s in the Bose glass 1T6| . The 
independence from phonons should also be reflected by the pre-exponential factor (Tq. In quasi-2d samples 
we expect it to be a number of order e^//i. This is conjectured from the analogy with the conduction at the 
mobility edge in (non-interacting) 2d systems undergoing a metal-insulator transition. For those systems 
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Fig. 2 Spectrum of the Bose glass. Due to the presence of Goldstone modes the superconducting phase has a continu- 
ous spectrum at all energies and temperatures, (a) At T = the spectrum in the Bose glass is discrete at low energies, 
reflecting the localization of bosons. The continuum starts at an non-extensive energy Ec for intermediate disorder 
Qc < g < g*. Above the mobility edge is extensive, (b) At finite temperatures, the discrete spectrum only survives 
in the regime with an extensive mobility edge. The continuum extends down to a; = everywhere below g^, due to 
inelastic scattering from thermal excitations in the continuum. 

the scaling theory of localization predicts a very weak temperature dependence of aoh/e'^ « 0(1). Note 
that such a prefactor is also natural to expect if the critical conductivity is to be recovered in the limit 
g Qc- Empirically one indeed consistently finds prefactors of this order 

The T dependence of the conductivity is entirely different for stronger disorder, g > g*. Since the 
mobility edge is at an extensive energy, the system will remain localized at low temperature. However, as 
T increases the number of inelastic scattering channels increases, and in some (possibly narrow) regime of 
disorder < g < goo one expects diffusion to be possible above a finite localization temperature li24ll2Tl 
T\oc{g) > 0. Approaching the latter from above, one expects a divergence of the effective activation energy, 
^Aig > gc,T T\oc{g)) oo. Finally, if the band width of boson states is finite and interactions are 
short ranged, there is a maximal possible inelastic scattering rate. When it falls short of the relevant level 
spacing in strong enough disorder, g > goo, the Bose glass will be localized at any temperature. 

We emphasize that the non-extensive mobility edge Ec is only well defined at T = 0. At finite T in- 
elastic collisions introduce a finite life time for excitations below Ec- According to the Thouless argument, 
diffusion is then already possible at energies E < Ec where the relevant level spacing is smaller than the 
inverse life time. This leads to an effective lowering of the "diffusion edge" T) < Ec{g, 0) with 

growing temperature. Such a tendency to superactivation was actually observed in recent experiments ||251 . 

Note that such a temperature dependence is also required to match the two regimes g < g^, and g > g^,. 
At fixed temperature, Ea increases monotonously with g, and diverges only beyond g». In particular 
EA{g*,T > 0) is finite, and it must be bigger than EA{g < g*,T). This is only consistent with EA{g 
g*, T = 0) = oo, if for fixed g close to the activation energy Ea decreases with increasing temperature. 

3 Discussion 

The above scenario of the Bose glass should be experimentally relevant at temperatures where weakly 
phonon-assisted variable range hopping can be neglected. In order to justify the bosonic description the 
parameters must be such that Ea is smaller than the depairing energy of electron pairs, restricting the 
validity of the picture to moderate magnetic fields. The experimentally most relevant part of the phase 
diagram is the conducting part of the Bose glass which exhibits purely electronic activation of Arrhenius 
type. It would be interesting to test this scenario by confirming that transport in that regime is indeed by 
bosons, and by complementing transport studies with probes which are sensitive to the spectrum, such as 
absorption measurements or low frequency AC conductivity. 
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One may worry that the notion of a mobility edge Ec makes an unreasonable assumption about the 
homogeneity of the disorder distribution tacitly implied when writing down the Hamiltonian ([T]i. In this 
context it is worthwhile noticing that the typical scale of Ec is of order t, and thus of order the bulk 
Tc of the material. The bosonic mobility edge thus lives on a much lower scale than mobility edges 
in fermionic Anderson problems, which scale with the Fermi energy. Accordingly, the sensitivity with 
respect to inhomogeneity is much smaller at the same temperatures. Except very close to the SIT, we thus 
expect the bottleneck of conduction to be due to the exponentially small probability of finding excitations 
in the continuum, even though the spectral weight at the mobility edge of the continuum will have certain 
spatial variations. 

In the discussion of simply activated transport channels we have assumed that the excitations at the 
mobility edge could carry charge. This seems plausible for the Bose glass, but is not necessarily true in 
more general systems close to a quantum delocalization transition. In Ref. 1261 a scenario for delocalized 
collective electronic vibrations in quantum electron glasses was considered. Those only transport energy, 
but no charge. If the low energy part of the continuum is entirely formed by such excitations one still 
expects purely electronic simple activation at the lowest temperatures ll26l . However, there will also be 
a significant range of higher temperatures where the optimal conduction channel is the variable range 
hopping of charged carriers, assisted by the delocalized vibrations (plasmons) of the electron glass. We 
suspect that a similar scenario applies to the high field regime of dirty superconducting films where most 
bosons depair into electrons, while the purely electronic nature of transport persists. 

In the present paper we emphasized the physical picture and the ideas leading to what we believe to be a 
rather generic phase diagram of localization phenomena close to a transition from a disordered insulator to 
a conducting phase. Calculations on the Bethe lattice, which are well controlled at high connectivity 1271 . 
confirm most qualitative properties of the phase diagram. Analytical tools allowing for a verification of the 
discussed scenario are presently being developed. 
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